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ABSTRACT 

This  paper  examines  the  crushing  charact eristics  of  a  web 
girder  subjected  to  a  local  in-plane  crushing  load.  First,  a 
simplified  model  of  the  post-buckling  deformation  of  the  web 
girder  is  developed.  Then  this  model  is  used  to  formulate  an 
upper-bound  solution  for  the  plastic  behavior  of  the  deforming 
web  girder.  Finally  small-scale  crushing  tests  of  extruded 
aluminum  double-hollow  forms  are  conducted  to  verify  the 
deformation  mode  and  the  upper-bound  solution.  The  web  girders 
of  interest  for  this  research  are  the  longitudinal  web  girders 
in  a  unidirectionally  stiffened  double  hull  (USDH)  structure. 
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1 .   INTRODUCTION 

The  ability  of  ships  to  withstand  damage- -whether  from 
collision,  grounding,  or  explosion--is  an  active  concern  of 
naval  architects  and  ship  owners.  In  1985,  the  Hitachi  Zosen 
Corporation  of  Japan  presented  the  concept  of  an  advanced  double 
hull  design  which  utilizes  longitudinal  girders  connecting  inner 
and  outer  hulls  to  form  a  cellular  structure  that  wraps  around 
the  girth  of  a  hull.111  Among  other  advantages,  this 
unidirectionally  stiffened  double  hull  (USDH)  structure,  shown 
in  Figure  1,  offers  the  potential  of  increased  protection 
against  local  damage  resulting  from  accidental  loading.  This 
potential  has  contributed  to  recent  interest  in  the  USDH 
structure  for  both  commercial  and  naval  ship  designs.121 

Damage  to  ships  involves  a  complex  interaction  of 
structural  members.  While  several  authors  have  developed  models 
of  ship  collisions,  research  concerning  ship  groundings  is  in 
the  initial  stages .  Most  ship  damage  models  have  dealt  with 
conventionally  framed  ships  which  utilize  a  gridwork  system  of 
large  transverse  web  frames  intersected  by  smaller  longitudinal 
stiffeners.  In  general,  for  both  collisions  and  groundings, 
these  transverse  web  frames  are  subjected  to  large  compressive 
in-plane  loads . 

McDermott  et  al.[3]  developed  a  mathematical  model  of  the 
various  stages  involved  in  a  minor  collision  of  longitudinally 
framed  oil  tankers.   This  model  was  applied  to  both  single  hull 


and  double  hull  tankers.  McDermott  assumed  that  the  resistance 
offered  by  a  buckled  web  frame  remained  constant  during 
deformation  of  the  web  frame  in  the  transverse  direction. 
Changt4]  in  his  discussion  of  McDermott' s  paper  recommended  that 
the  post-buckling  resistance  of  web  frames  be  further 
investigated  because  many  structural  members  offer  decreased 
resistance  after  an  ultimate  strength  has  been  reached. 

Wierzbicki  et  al . [5]  investigated  the  extent  of  damages 
involved  in  high  energy  groundings  of  tankers  for  a  study 
conducted  for  the  Joint  M. I . T . -Industry  Program  on  Safe  Tankers. 
For  this  study,  a  simplified  model  was  developed  to  predict  the 
energy  dissipated  by  the  crushing  of  transverse  web  frames. 

This  model  utilized  the  assumption  that  the  deformation  of 
the  transverse  web  frame  and  the  hull  plating  which  acts  as  a 
flange  could  be  separated  into  two  mechanisms.  It  was  assumed 
that  the  web  detaches  from  the  flange  and  rolls  into  two 
cylinders  of  equal  radii,  one  on  either  side  of  the  point  of 
load  application.  It  was  assumed  that  the  hull  plating  deforms 
into  a  trough  shape.  A  sketch  of  this  model  is  shown  in  Figure 
2.  Geometric  compatibility  of  the  connection  between  the  web 
and  the  attached  plating,  however,  was  not  maintained. 

Instead  of  simplifying  the  post-failure  resistance  of 
deforming  webs  by  modelling  the  resistance  as  a  constant,  the 
plastic  behavior  of  the  web  may  be  described  by  a  non-linear 
spring  in  a  global  model  of  a  ship  structure . 


For  example,  the  USDH  structure  can  be  considered  to  be  a 
sandwich  structure  with  a  definable  crushing  strength  of  the 
core,  a  finite  shear  strength  in  the  direction  of  the 
longitudinal  web  girders,  and  negligible  shear  strength  in  the 
transverse  direction  normal  to  the  web  girders .  This  sandwich 
structure  can  be  modelled  as  rigid-plastic  membranes 
(representing  the  inner  and  outer  hull  plating)  connected  by 
vertical  (crushing  strength)  and  diagonal  (shear  strength)  non- 
linear springs . [S]   This  model  is  shown  in  Figure  3 . 

The  objective  of  this  study  is  to  develop  the  load- 
deformation  characteristics  for  longitudinal  web  girders 
subjected  to  local  in-plane  crushing  loads.  These 
characteristics  can  then  be  used  to  describe  the  vertical 
springs  of  the  model  described  above . 

The  crushing  of  longitudinal  web  girders  presents  a  problem 
similar  to  that  of  crushing  tubular  members.  Several  studies 
have  been  conducted  which  have  explored  the  large  deflection 
plastic  behavior  of  thin  walled  tubular  members.  Reference  [7] 
contains  a  review  of  several  of  these  studies .  Much  of  the 
analyses  involving  crushing  of  tubular  members  is  focused  on 
crushing  due  to  axial  loads.  This  class  of  problems  lends 
itself  to  simple  modelling  due  to  the  axi-symmetrical  cross- 
section  of  tubular  members . 

The  crushing  of  longitudinal  girders  in  the  USDH  structure, 
however,  is  similar  to  crushing  of  tubes  subjected  to  lateral 
loads  rather  than  axial  loads.   De  Oliveira  et .  al .  [81   present 


an  analysis  of  crushing  of  circular  tubes  due  to  lateral  loads. 
This  class  of  problems  is  more  difficult  to  analyze  due  to  a 
lack  of  symmetry  and  due  to  the  fact  that  a  tube's  cross- 
sectional  shape  changes  during  deformation.  If  the  crushing 
load  is  localized  rather  than  distributed,  the  distortion  of  the 
cross-sectional  shape  occurs  locally  and  propagates.  This 
propagation  of  damage  is  difficult  to  predict  with  out 
experimental  results . 

The  problem  under  consideration  for  this  study  resembles 
crushing  of  a  multiple  box  column-type  structure  subjected  to  a 
localized  lateral  load.  Although  de  Oliveira's  work  focused  on 
crushing  of  circular  tubes,  the  same  approach  of  using  simple 
methods  of  structural  plasticity  applied  to  a  simplified 
deformation  model  is  used  for  this  study. 

First,  models  for  two  modes  of  deformation  of  a  web  girder- 
-local  denting  and  global  bending--are  developed  based  on  some 
photographic  evidence  of  the  deformation  of  fully  clamped  plates 
that  were  subjected  to  in-plane  impact  loading.  Next,  upper- 
bound  solutions  for  the  plastic  behavior  of  deforming  web 
girders  are  formulated  for  each  mode.  Finally,  theoretical 
predictions  based  on  the  upper-bound  solutions  are  compared  to 
the  results  from  a  limited  number  of  small-scale  crushing 
experiments . 


2.   PROBLEM  DEFINITION 

2.1.  UPPER-BOUND  THEOREM:  An  upper-bound  approach  will  be 
used  to  find  the  crushing  characteristics  of  a  web  girder 
subjected  to  a  local  in-plane  crushing  load.  Calladine  [91 
provides  a  concise  statement  of  the  Upper-bound  Theorem: 

"If  an  estimate  of  the  plastic  collapse  load  of  a  body 
is  made  by  equating  internal  rate  of  dissipation  of 
energy  to  the  rate  at  which  external  forces  do  work  in 
any  postulated  mechanism  of  deformation  of  the  body, 
the  estimate  will  be  either  high,  or  correct." 

The  postulated  mechanisms  of  deformation  to  which  Calladine 
refers  to  must  be  kinematically  admissible. 

Two  modes  of  deformation  are  postulated  for  this  study. 
The  first  mode  involves  local  indentation  of  the  web  girder. 
The  second  mode  involves  the  development  of  a  three  hinge 
mechanism  to  induce  global  bending  of  the  web  girder  between 
transverse  bulkheads.  These  postulated  deformation  mechanisms 
are  described  in  sections  3  and  4 .  Assumptions  concerning  the 
geometry  and  boundary  conditions  of  the  deforming  body,  the 
material  properties  of  the  body  and  the  applied  external  loads 
for  this  solution  are  described  in  the  following  sections. 


2.2.  STRUCTURAL  UNIT:  Figure  4  shows  a  section  of  the  USDH 
structure  consisting  of  rectangular  cells.  This  cellular 
structure  runs  between  transverse  bulkheads  and  wraps  around  the 
girth  of  a  hull.  Due  to  symmetry  it  is  sufficient  to  define  a 
structural  unit  which  consists  of  one  longitudinal  web  girder 
with  the  attached  inner  and  outer  hull  plating.  This  structural 
unit  is  shown  in  Figure  5. 

The  depth  of  the  web  girder  is  "H"  and  the  thickness  of  the 
web  girder  is  "tw"  .  Because  the  inner  and  outer  hull  plating 
act  as  fully  effective  flanges,  the  breadth  of  the  attached 
plating  is  equal  to  the  girder  spacing,  "B" .  The  thickness  of 
the  hull  plating  is  "tp"  .  The  length  of  the  web  girder  is 
determined  by  the  distance  between  transverse  bulkheads.  This 
length  is  much  greater  than  the  girder  depth  and  spacing.  The 
structural  unit,  therefore,  can  be  assumed  to  be  effectively 
infinite  in  the  longitudinal  direction. 

Due  to  symmetry,  the  attached  hull  plating  flanges  are 
constrained  from  in-plane  movement  in  the  X-direction  along  the 
longitudinal  edges  (X  =  ±  B/2)  .  Both  ends  of  the  structural 
unit  are  considered  to  be  fully  clamped  at  the  transverse 
bulkheads  (Y  =  ±  L/2) . 

2.3.  LOADING  CONDITION:  The  web  girder  is  considered  to  be 
subjected  to  a  local  crushing  load  such  as  that  delivered  by  a 
knife-edge  indenter  or  a  rigid  punch  normal  to  the  longitudinal 
axis  of  symmetry.    This  loading  is  quasi-static,  therefore 
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dynamic  effects  are  ignored.   The  applied  load  has  a  magnitude 
of  PB.   This  load  causes  the  flange  to  deflect  a  distance  A 
relative  to  the  un-deformed  region  of  the  flange.   The  rate  of 
deflection  is  denoted  by  A. 


2.4.  MATERIAL  PROPERTIES  AND  ASSUMPTIONS:  The  material  is 
assumed  to  be  rigid/plastic  isotropic  time -independent  material 
which  can  be  described  by  the  flow  stress  <J0.  Strain  hardening 
effects  may  be  included  by  adjusting  the  value  of  the  flow 
stress  based  on  the  average  plastic  strains  of  the  deforming 
regions  during  the  deformation  process.  The  average  plastic 
strain  may  be  obtained  by  integrating  the  strain  over  the  length 
of  the  deforming  region: 


=  -  f  e(z)  dz      .  (1) 


€avs 

0 


Bending  stresses  and  membrane  stresses  are  modelled  as 
uncoupled.  Deforming  regions  which  experience  high  bending 
stresses  are  considered  to  have  negligible  membrane  stresses. 
Likewise,  regions  which  develop  significant  membrane  stresses 
are  assumed  to  have  negligible  bending  stresses . 
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2.5.  LIST  OF  NOMENCLATURE:   This  section  provides  a  list  of 
the  symbols  and  parameters  used  in  this  paper. 


HULL  PARAMETERS 

H: 
B: 
t„: 

V 

L: 

H*: 

B*: 

f: 

L*:         L  /  H; 


Depth  of  web; 

Longitudinal  girder  spacing; 

Thickness  of  web; 

Thickness  of  inner  and  outer  hull  plating 

Distance  between  transverse  bulkheads; 

H  /  tw; 

B  /  tp; 

tp  /  t„; 


MATERIAL  PARAMETERS 


Eult 

M0 
M0 
Mu 
Ne 
NL 


Flow  stress; 

Ultimate  tensile  yield  stress; 
Ultimate  tensile  strain; 

Fully  plastic  bending  moment  for  the  web; 
Fully  plastic  bending  moment  for  the  flanges; 
(Ou  tK2)/4 

Fully  plastic  membrane  force  for  the  web; 
Fully  plastic  membrane  force  for  the  flanges; 


GEOMETRIC  PARAMETERS 


A: 

A*: 
5: 

C: 

a: 


Crushing  depth; 

A  /  H; 

Local  crushing  depth  for  global  bending  mode; 

Vertical  extent  of  damage  for  local  denting  mode; 

5  /  H; 

Angle  between  first  and  second  hinge  lines; 

Angle  between  second  and  third  hinge  lines; 
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GEOMETRIC  PARAMETERS,  CONTINUED 


e0: 

♦ 
Y 

uf: 

e£: 
ex: 


Angle  of  indentation  due  to  crushing; 

Angle  between  regions  connected  to  the  i'th  hinge 

line; 

Length  of  i'th  hinge  line; 

Angle  between  the  first  and  third  hinge  lines; 

Angle  between  center  line  and  third  hinge  line; 

Angle  of  depression  in  global  bending  mode; 

Stretching  of  flange; 

Maximum  stretching  of  web; 

Average  strain  for  flange; 

Maximum  strain  for  web; 

Average  strain  for  structural  unit; 

Crushing  force . 
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3.   LOCAL  DENTING  MODE 

3.1.  SIMPLIFIED  MODEL:  Figure  6  shows  a  photograph  of  a 
fully  clamped  plate  which  was  deformed  by  a  crushing  load 
applied  by  a  rigid  punch  to  the  upper  edge  of  the  plate.  This 
photograph  is  from  full  scale  experiments  conducted  by  S. 
Shimizu  of  Shimushu  University  in  Japan .  [101  Based  on  this 
photograph  and  other  photographs  of  damage  incurred  by  tankers 
during  grounding  accidents,  the  model  for  the  local  deformation 
mechanism  shown  in  Figure  7  was  developed. 

The  deformation  of  the  web  involves  both  bending  and 
stretching.  The  bending  of  the  web  is  caused  by  the  formation 
of  plastic  hinge  lines  which  are  shown  in  Figure  7.  The  web 
deforms  symmetrically  about  the  point  of  application  of  the 
load.  Six  stationary  plastic  hinge  lines  form.  Due  to 
symmetry,  only  the  three  hinge  lines  on  one  side  of  the  vertical 
plane  of  symmetry  will  be  described.  The  first  hinge  line  forms 
at  the  junction  between  the  web  and  the  upper  flange.  The 
second  hinge  line  forms  at  an  angle  a  below  the  first  hinge 
line.  The  third  hinge  line  forms  at  an  angle  |3  below  the  second 
hinge  line. 

The  photographs  examined  in  developing  this  model  showed 
the  length  of  the  deformed  region  in  the  longitudinal  direction 
to  be  approximately  twice  the  length  in  the  vertical  direction. 
This  characteristic  was  included  in  the  model  by  fixing  the 
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angle  between  the  upper  edge  of  the  undeformed  web  and  the  third 
hinge  line  at  tc/4.  The  vertical  length  of  the  deformed  zone  is 
taken  as  a  parameter  £. 

3.2.  DEFORMATION  SEQUENCE:  Figure  8  depicts  the  geometry 
of  the  deformed  web.  The  applied  force  PB  causes  the  upper  edge 
of  the  web  to  indent  a  distance  A.  The  angle  of  indentation  is 
0O.  As  A  increases  the  second  hinge  moves  out  of  plane.  The 
angle  between  the  first  and  third  hinges,  <X>,  decreases.  The 
angle  between  the  vertical  plane  of  symmetry  and  the  third  hinge 
line,  <J>,  remains  constant  at  7i/4. 

The  shaded  area  in  Figure  8  illustrates  the  stretching  in 
the  deforming  web.  The  attached  flange  deforms  with  the  web  as 
shown  in  Figure  9.  Again,  the  shaded  area  illustrates  the 
stretching  in  the  deforming  flange. 

3.3.  LOAD-DEFLECTION  RELATIONSHIP:  Following  the  Upper- 
Bound  Theorem,  the  load  deflection  relationship  is  derived  by 
equating  the  rate  of  external  work  performed  in  deforming  the 
web  girder  to  the  rate  of  internal  energy  dissipation.  In  the 
following  sections,  this  derivation  is  discussed.  First  the 
rate  of  internal  energy  dissipation  is  derived.  Next,  the  rate 
of  external  work  is  formulated.  Then,  these  two  rates  are 
equated  to  yield  the  load-deflection  relationship. 
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3.3.1.    RATE  OF  INTERNAL  ENERGY  DISSIPATION:    The 

total  rate  of  energy  dissipation  consists  of  two  components-- 
that  contributed  from  bending,  EB,  and  that  contributed  from 
stretching,  F^. 


ETOTAL  EM       +  EB 


(2) 


3.3.1.1.   BENDING  ENERGY  RATE:    In  general  the 

bending  energy  rate  term  consists  of  energy  dissipated  by  the 

continuous  deformation  field  and  energy  dissipated  by  the 
discontinuous  deformation  field. 

EB     =     f      M0k     dS      +       t     M0   [i±]    ld  (3) 

•*  s  1=1 


where  S  is  the  surface  area  of  the  deforming  region;   K  is  the 
rate  of  change  of  curvature;  Qt   is  the  rate  of  rotation  in  the 
i'th  hinge  line;  and  1L   is  the  length  of  the  i'th  hinge  line. 

For  the  local  denting  model,  however,  the  deformation  is 
simplified  by  modelling  all  plastic  hinges  as  stationary  hinges. 
There  is  no  continuous  deformation  field.   The  bending  energy 
rate  term  simplifies  to: 

EB    =  2  M0    {    [6X]  1,  +   [4a]  I2  +   [43]  I,  }   ,        (4) 
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where  M0  is  the  fully  plastic  bending  moment: 


m     a^_ti  (5) 


The  assumption  is  made  that  all  hinge  lengths  remain  constant 
during  the  deformation  process .  The  hinge  lengths  can  be 
determined  from  the  geometry  of  the  model  shown  in  Figure  6 . 


A  =  C 

12     =  C  /  cos (a)  (6) 

i3  =  ^  C 


The  rotation  rates,  Q^,    for  each  hinge  line  can  be  related 
to  the  angle  <X>  between  the  first  and  third  hinge  lines.   For 

each  hinge  line,  a  section  normal  to  the  hinge  line  is  used  to 
determine  the  rotation  rate. 

From  the  geometry  of  the  model,  the  following  can  be 
determined: 

~~„  a            cos  B  -  cos  a   cos  $  ni 

COS  0,    =   H ; -j- ,  I  ') 

sin  a  sin  v 


cos  6,  =  cos  a  -  cos  p  cos  fl  (8) 

sin  p  sin  <D 
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cos  0,  =   -  cos  a   cos  g  +  cos  *  (9) 

sin  a  sin  p 


Using  the  relationship  that: 


-£     [arccos(u)]     =  u  (10) 

dt  /(l  -  u2) 


with    (7) ,     (8) ,    and    (9) ,    the   rotation   rates    can  be   shown   to  be 
A      =  -  *(coscc  -  cosp  cos*) 


sin*  sjC1  +  Cz  cos*  -  cos2* 


(11) 


a  4>    sin*  ,„_. 

92  — ==========       '  (12) 

^Cx  +  C2  cos*  -  cos2* 


A      _         -  *     (cosP  -  cosa  cos*  ) 


sin*    Jc^  +  Cz  cos*  -  cos2* 


(13) 


where 


Cx    =     (sin2a  -  cos2p) 
C2     =     (2  cosa  cosp) 


(14) 


By  noting  that  the  sum  of  the  depression  angle,  0O,  and  <I> 
remains  constant  and  equal  to  rc/4,  COS  (O)  ,  SIN(O),  and  O  can  be 
related  to  A  and  A. 
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cos<D     =     H     Sj  +  A)  (15) 

2     v/C2  +  A2 


sin*     =     V?     (C  "  A)  (16) 

2      7C2  +  A2 


&   =   — zA_L_  (iv) 

(C2  +  A2) 


Substituting    (11),     (12),     (13),    and    (6)     into    (4)    yields 


E  --7.1M  j»  (cos2a-cosacospcos*)  +v/2"(costtcosP-cos2acos*)  +2sin2* 

(18) 


<.B--Z(M0*, 

coso  sin*  yjC^  +  C2cos*  -  cos2* 


-2  C  MQ  * 


(C3  -  C4  cos*   -  cos2*) 


coso  sin*  SJC1  +  C2  cos*  -  cos2* 


(19) 


where 


C3     =     (cos2a  +  v/2cosacosP   +  1) 

C4     =     (N/2cos2a  +  cosacosp)        .  (2°) 
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3.3.1.2.  MEMBRANE  ENERGY  RATE:  Now  consider  the 
membrane  energy  rate  term  of  (2)  .  Both  the  attached  upper 
plating  (or  flange)  and  the  web  stretch  while  deforming. 


Therefore 


*         ♦  t  ... . 

EM  total      ~      EM  flange      +      EM  web        ■  l*1/ 


The  rate  of  membrane  energy  dissipation  is  defined  as: 

EM    =  f     N0k  dS    =  f     N0k  dr)dl        ,  (22) 

J  s  J  s 


where  e  is  the  velocity  strain  in  the  deformed  co-ordinates 
(£,T|)  and  N0  is  the  fully  plastic  membrane  force,  defined  below: 

4  =  ^f  •  <23> 

dl 


N0    =  o0t       .  (24) 

Integration  of  (22)  yields: 

EM    =  2N0    /"iSdr,   ,  (25) 
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where  2u  is  the  velocity  with  which  material  points  on  either 
side  of  the  centerline  of  the  shaded  regions  shown  in  Figures  8 
and  9  move  apart  relative  to  each  other. 

First,  consider  the  stretching  involved  in  deforming  the 
flange.  Figure  9  shows  the  trough  shape  of  deformation  that  the 
flange  assumes.  Although  the  stretching  of  the  flange  occurs 
across  the  entire  deformed  region,  the  shaded  area  in  Figure  9 
is  used  to  approximate  the  average  strain  involved.  The  amount 
of  stretching  is  constant  in  the  ^-direction,  normal  to  the 
longitudinal  girder. 

From  the  geometry  of  the  model,  it  follows  that: 

u   =     £    -  r     =  V^  C _  r       t  (26) 

f     cos0o        (cos*  +  sin*) 


ut 


■y/2   *  C  (cos*  -  sin*)  (27) 

(1+2  cos*  sin*) 


Therefore, 


EM  flange      =      2    N0    f     U  dl\ 


(28) 


=  -2^2  W0  B  C  i     /cos*   -sin*) 
°  (1  +  2cos*sin*) 


21 


where  N0  is  the  fully  plastic  membrane  force  for  the  flange: 


No    =  o0  tp 


(29) 


Now,  consider  the  stretching  involved  in  deforming  the  web. 
Again,  the  shaded  region  in  Figure  8  can  be  used  to  approximate 
the  average  strain  for  the  web.  The  amount  of  stretching  varies 
along  the  depth  of  the  web.  The  displacement  u£  at  the  web 
flange  connection  and  the  maximum  displacement  ux  at  the 
intersection  of  the  second  hinge  lines  are  used  to  determine  the 
average  displacement  as  shown  in  Figure  10. 


u 


web.avg 


/ 


Ctano 


/ 


c 

Ctano 


Uf   + 


(ux  -  uf)  z 
Ctancc 


ux  (C  -  z) 
C(l  -  tana) 


dz 


dz 


(30) 


Uwebavg     =      \    ux     +      \    tano    uf 


(31) 


web  avg 


2       x 


-1 

—  tana  u< 


(32) 
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From  the  geometry  of  the  model,  the  following  can  be  determined 

x  [  (y/2  cosa  -  cosP)sin*  +  cosP  cos*  -  cosaj    /33V 

x  [cosP  (sin*  -  cost)  +  cosa] 


u      -     -x/2  *  C    (cosacosp  -cos2acos*)  ^34j 

[cosp  (sin*  -cos*)  +  cosa]2 


Substituting  (27)  and  (33)  into  (31)  and  then  substituting 
(31)  into  (25)  yields: 


JM  web 


=     -t/2  *  C2  N      (  (cosacosp  -  cos2acos*) 

[cosP(sin*  -  cos*)    +  cosa]2        .3_. 

sina(cos*  -  sin*) 


cosa    (1  +2  cos*  sin*) 


) 


3.3.2.   RATE  OF  EXTERNAL  WORK:   The  rate  of  work  done 
in  deforming  the  structural  unit  is : 

*«.  -  p*i    •  <36> 

Rearranging  (17)  and  substituting  this  into  (35)  gives: 


»m  -   -  pb  *  JSL^L     .  <37, 
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3.3.3.  PHYSICAL  LIMIT:  For  the  local  denting  mode, 
a  is  modelled  as  a  parameter.  The  choice  of  a  affects  not  only 
the  magnitude  of  the  crushing  force  for  a  specified  crushing 
depth,  but  also  the  range  of  freedom  for  the  model.  The  degree 
of  freedom  for  the  denting  mode  is  exhausted  when  the  upper 
flange  makes  contact  with  the  "upper  flap"  of  the  web,  the 
region  between  the  first  and  second  hinge  lines.  When  this 
occurs,  0!  =  7t/2  and  : 


'  A*\  _  cos2a  +  2  cosP  y/cos2g  -  cos2P 


(38) 
I  C  )limit  2   cos2p  -  cos2a 


For  approximately  a  <,  n/12,  the  linear  assumption  discussed 
in  Appendix  A  is  valid.  The  load-deflection  relationship  in 
dimensionless  terms  is: 


PB  2[{KX   -  K2)    -    (K2   +  2)i|r] 


M0 


COSO(l  -  !|f)  J(K3    -    KJ     +  (tf4  -  2)!|l 

2^/2  CH*    [(K5   -   #6)  -  K^]  (39) 


cos2a  -  K.ty 


4 

2' 


+  i|i  [4  C*#*  tana  +  8  B*t* 


where  the  dimensionless  parameters  are: 


A*  =  A   ;     {*     =  I   ;     H*     =  JL      ; 

H  H  t. 


vs 


B*  =  -B.   ;  t*  =  ie 


24 


(40) 


and 


A* 
i*r  =  (41) 


The  constants  Kx  through  K7  are  functions  of  a  and  (3  which  are 
listed  in  Appendix  A  (A6) . 


3.3.4.  INITIAL  PHASE  ENVELOPE:  For  the  local  denting 
mode,  it  is  postulated  that  the  longitudinal  and  vertical  extent 
of  deformation,  £,  is  initially  zero  but  spreads  during  the 
deformation  process  until  the  deformed  region  extends  the  depth 
of  the  web  {C,   =   H)  . 

Figures  11  and  12  show  the  local  denting  solution  for 
values  of  C,*  increasing  from  C,*  =  0.1  to  C,*  =  1.0.  The  parameters 
for  the  hull  geometry  (H*  =  100;  B*  =  75;  t*  =  1.0)  in  both 
Figures  11  and  12  are  approximate  values  for  a  40,000  DWT  USDH 
tanker.  Figure  11  depicts  the  load-deflection  curve  for  a  = 
tc/12.  The  load-deflection  curve  for  a  =  7i/8  is  shown  in  Figure 
12. 

Both  Figures  show  envelopes  for  the  initial  phase  of 
deformation.  Due  to  the  complexity  of  the  load-deflection 
equation  (39)  it  is  difficult  to  analytically  determine  this 
envelope.  The  initial  phase,  however,  may  be  described 
empirically  by  a  function  of  the  form: 
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£  -  c°{Ay 


(42) 


where  the  constant  C0  and  the  exponent  "m"  depend  on  both  the 
value  of  alpha  and  the  geometric  parameters  describing  the 
structure. 

3.3.5.  STRAIN  AND  STRAIN  HARDENING:  From  the 
geometry  of  the  model  the  strain  averaged  over  the  entire 
deforming  area  can  be  shown  to  be : 


'&vg 


B't*   +  -£lr 
H* 


'  +  -j—  tana 
2h* 


B*t 


€f  + 


2h* 


.       (43) 
) 


where 


_      y/2   -  cos$   -  sinfl 
(cos$  +  sinfl>) 


(44) 


and 


(v/2cosa  -  cosP)sin$  +  cosPcosfl  -  cosce 
[cosP(sin$  -  cos*)    +  coso] 


(45) 


Strain  hardening  effects  can  be  accounted  for  by  adjusting 
the  value  of  a0  based  on  the  average  plastic  strain.  For 
example,  if  the  stress-strain  relationship  for  a  given  material 
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can  be  described  by  a  power-law  relationship  similar  to  the 
Ramberg-Osgood  formula  [11)  : 


°o   =  °ult 


t €    \n 


V  ^ult, 


(46) 


then,  the  following  refinement  of  the  denting  model  could 
be  made : 


/  e         \ 
-avg 


M, 


ult 


\£ultj 


a     p 


(47) 


where 


M      = 


°ult    tw2 


(48) 


and  PB/Mo   is   defined  by    (39)    or    (A5) . 
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4.   GLOBAL  BENDING  MODE 

4.1.  SIMPLIFIED  MODEL:  This  second  mode  of  failure 
involves  generation  of  a  three  hinge  beam-bending  mechanism, 
shown  in  Figure  13.  For  simplicity,  the  crushing  load  is 
applied  at  mid- span  and  the  length  between  "hinge  points"  is 
equal.  This  mechanism  could  be  adjusted  to  account  for  unequal 
lengths  between  hinge  points . 

The  "hinge  points"  are  modeled  as  shown  in  Figure  14.  This 
model  is  similar  to  the  local  indentation  model  except  there  is 
no  extension  of  the  flange.  The  vertical  extent  of  deformation, 
£,  extends  the  entire  depth  of  the  web;  therefore,  C,  =  H.  Also, 
a  fourth  hinge  line  forms  along  the  bottom  flange  at  the 
centerline  of  the  "hinge  point"  mechanism. 

4.2.  DEFORMATION  SEQUENCE:  Figures  13  &  14  show  the  web 
deforming  in  the  global  mode.  The  applied  force  PB  causes  the 
legs  of  the  three  hinge  mechanism  to  rotate  an  angle  y.  This 
causes  the  angle  <j>  for  each  "hinge  point"  mechanism  to  decrease. 
The  second  hinge  line  moves  out  of  plane  and  the  angle  between 
the  first  and  third  hinge  lines,  O,  decreases.  This  results  in 
a  local  deflection,  5,  at  the  central  "hinge  point"  mechanism. 
This  local  5  is  part  of  the  total  deflection  A  as  shown  in 
Figure  13.  The  shaded  area  in  Figure  15  illustrates  the 
stretching  of  the  web  involved  in  deforming  the  "hinge  point" 
mechanism. 
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The  range  of  A*  is  limited  to  A*  <,  0.5.  For  a  ship  which 
may  have  a  web  girder  depth  on  the  order  of  2  meters,  A*  =  0.5 
is  a  significant  displacement  of  the  hull.  Because  the  range  of 
A*  is  limited,  the  extension  of  the  two  "legs"  of  the  global 
mechanism  is  neglected. 

4.3.  LOAD-DEFLECTION  RELATIONSHIP:  The  application  of  PB 
causes  the  three  "hinge  point"  mechanisms  to  form.  The  "hinge 
point"  mechanisms  are  separated  by  a  distance  (L/2  -  H)  .  The 
central  "hinge  point"  is  designated  the  primary  mechanism.  The 
two  outer  "hinge  points"  are  designated  secondary  mechanisms. 
The  subscripts  "p"  and  "s"  are  used  to  designate  the  primary  and 
secondary  mechanisms,  respectively.  As  the  central  mechanism 
deflects,  the  two  "legs"  of  the  overall  structure  rotate  an 
angle  Y-  From  Figures  13  and  14,  the  angles  <j>p  and  <|>.  can  be 
related  to  y. 


4>p  =  |  -  Y  (49) 


and 


♦.  -  i  -  i    •  <50) 


29 


In  the  following  sections,  the  internal  rate  of  energy 
dissipation  for  the  primary  mechanism  is  formulated  as  a 
function  of  y  and  y,  fp(y,y)  .  In  the  global  bending  mode,  energy 
is  dissipated  not  only  by  the  primary  "hinge  point"  mechanism 
but  also  by  the  two  secondary  "hinge  point"  mechanisms.  The 
total  energy  dissipated,  therefore,  is: 


Etotai    =    fpiy.i)  +  2fs(Y,Y) 


(51) 


where  fp  and  f,  are  the  functions  formulated  in  the  following 
sections  for  the  total  rate  of  internal  energy  dissipation  for 
the   primary   and   secondary    "hinge   point"   mechanisms, 
respectively.   The  functions  fp  and  f.  are  related  by  (4  9)  and 
(50)  so  that: 

'.<Y.f>  ■  f,(i.$\       •  <52> 


4.3.1.    RATE  OF  INTERNAL  ENERGY  DISSIPATION:    The 

deformation  of  the  web  girder  in  the  global  bending  mode,  like 
that  in  the  local  denting  mode,  involves  both  bending  and 
stretching.  The  expression  for  the  internal  rate  of  energy 
dissipation,  therefore,  is  the  same  as  stated  by  (2) . 
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4.3.1.1.  BENDING  ENERGY  RATE:  Like  the  local 
denting  model,  the  bending  model  consists  of  six  stationary 
hinge  lines  in  the  web.  In  addition,  one  hinge  line  forms 
across  the  bottom  flange  as  shown  in  Figure  13.  There  is  no 
continuous  deformation  field.  The  bending  energy  rate  term, 
therefore,  is: 

EB    =  2M0  {  lx  [6lp]  +  12  [62p]  +  I,  [4     ]    }   +  W0  I4  [64p]  (53) 


where  M0  is  the  fully  plastic  bending  moment  for  the  flange: 

m  =  ^jL  (54) 


The  assumption  is  made  that  the  lengths  of  the  hinge  lines 
remain  constant.  From  the  geometry  of  the  model,  shown  in 
Figure  14,  the  hinge  line  lengths  can  be  determined: 


(55) 


h 

= 

H 

h 

= 

H/cosa 

*3 

= 

s/2  H 

1* 

= 

H(l    - 

tana) 

Is 

= 

B 

First,  consider  the  fourth  hinge  line.   The  rotation  rate 
of  the  fourth  hinge  line  is  simply: 

e4P  =  2  {  (56) 
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Now  consider  the  other  hinge  lines.  The  rotation  rates, 
[QL]  ,  for  the  first,  second  and  third  hinge  lines  are  related  to 
the  angle  O  between  the  first  and  third  hinge  lines  by  the  same 
relationships  used  for  the  local  denting  model-- (11) ,  (12)  and 
(13)  .  This  is  due  to  the  similar  geometry  of  the  two  models. 
The  relationship  of  O  to  A,  however,  differs  from  that  used  for 
the  local  denting  model . 

To  relate  Op  to  y,    the  following  geometric  relationships, 
which  are  evident  from  Figure  15,  are  used: 

*p  -  (f  "  *p)  "  flo  <57) 


cos0o    =    y/2  sin<J>_  (58) 


sin60    =    ^cos2<|>p  -  sin24>p  (59) 


cos$p  =  i—  {   (l-sin2y)    +   (cosY+sirvy)  v/sin2y  }  (6°) 


sin*p  =  &{  cos2Y  -   (cosY-sinY)/sTn2v  }  («D 


32 


4_= 


P 

-2Y{sin2y^sin2Y-sin22Y+cos2Y  (  — -sin2Y)) 

\/sin2Y-sin22Y  vl_sin2Y+2sin22Y+2cos2Y  \/sin2Y-sin22Y 
In  Appendix  B,  the  EB  term  is  presented  as  a  function  of  y. 

4.3.1.2.  MEMBRANE  ENERGY  RATE:  The  rate  of 
membrane  energy  rate  term  is  defined  by  (22) -(25)  .  For  the  web, 
the  maximum  displacement,  ux,  shown  in  Figure  15,  is  given  by 

(32)  with  £  =  H. 

The  average  displacement,  uwaTg  for  the  shaded  region  in 
Figure  15  is: 

1  (63) 

;;  s    77 

web.avg  j     "x 

Differentiation  of  (63)  and  substitution  of  this  into  (25) 
yields : 

♦      _  •      ,      coscccosB  -  cos2cccos4>  ,_.. 

E„    =  -y/2   $_  iV0  H2   - - S —   .   (64) 

p  \  cosP (sin*p  -  cos$p)  +  coso  )2 

4.3.1.3.  SECONDARY  MECHANISMS:  The  internal 
rate  of  energy  dissipation  for  the  secondary  mechanisms  follows 
the  same  function  as  that  for  the  primary  mechanism  evaluated  at 
y/2 .  Both  the  bending  energy  rate  term  and  the  membrane  energy 
rate  term  for  the  secondary  mechanisms  are  listed  in  Appendix  B. 
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4.3.2.     RATE  OF  EXTERNAL  WORK:    The     rate     of 
external  work  performed  in  deforming  the  web  girder  is : 

W    *     =  PD  A   .  (65) 

The  local  deflection,  5,  at  the  central  mechanism  can  be 
described  as  a  function  of  y  by  considering  the  geometry  of  the 
model . 


6  =  JH2  -   2if2sin2<|>p  =  H  v^sin2Y  <66) 


The  total  deflection,  A,  is  the  sum  of  the  local  deflection  and 
the  deflection  due  to  the  rotation  of  the  two  legs  of  the  entire 
mechanism. 


A  =  fi  +  IL   -   #j  tany  (67) 


Therefore 


A  =  6  +  y  (—   -  H\   sec2y 


(68) 


Differentiation  of  (66)  and  substitution  of  the  result  into 
(68)  gives: 

1  _  y  (  H(cos2Y-sin2y)  cos2y  +  (L  -  2H)   Vcosysiny  }     /ggx 
2  cos2y  v'cosysinY 
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Substitution  of  this  result  into  (65)  gives: 


*«t  ■  (PBHy)   ZT     ,  (70) 


where 


Zt  ~~  in 


IT  V    ,  (71) 

(cos2y-sin2y)  cos2y  +  —  -  2  Vcosysiny 
2  cos2y  v/cosysiny 


4.3.3.  PHYSICAL  LIMIT:  As  in  the  local  denting 
mode,  a  is  modelled  as  a  parameter  for  the  global  bending  mode 
also.  The  choice  of  a  affects  not  only  the  magnitude  of  the 
crushing  force  for  a  specified  crushing  depth,  but  also  the 
range  of  freedom  for  the  model .  The  degree  of  freedom  for  the 
bending  mode  is  exhausted  when  the  connection  of  the  web  to  the 
flange  moves  out  of  the  Z-Y  plane.   This  occurs  when: 


V  limit    =  -T    arctan 


4  ^   sin2p 


/  2   \1 

1_  sin2a  2  (72) 
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4.3.4.     GLOBAL  LOAD -DEFLECT ION  RELATIONSHIP:     In 

Appendix  B,  the  external  rate  of  work  is  equated  to  the  total 
rate  of  internal  energy  dissipation  for  both  the  primary 
mechanism  and  the  two  secondary  mechanisms.  The  dimensionless 
equation  for  PB/Mo  is: 


B  _ 


M„ 


4Z. 


C3-C4cos$p-cos2G>p 


2J2H*  (cosocosB-cos2acos4t 

+ ! 2 

cosasinGp^q+qcosflp-cos2^      (cosP  (sin*p-cos*p)  +cosa) 


8Z. 


C3-C4cos*s-cos2$s 


2y/2H*  (cosacosp-cos2acos*s) 
■  cosasinfc^q  +  qcos^-cos2^       (cosP  (sin*s-cos*s)  +cosa) 2 


4B*t*3 
ZjH* 


where 


(73) 


A*     =    /2    (Y)  2     +      -=r   -1 


L* 


(74) 


and  the  a  new  dimensionless  parameter  is  introduced: 


l-  =  k 

H 


(75) 
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The  expression  for  ZT  is  given  in  (71)  .  The  expressions  for 
Zp,  COS  <X>p,  and  SIN  <E>p  are  given  by  (B2)  ,  (B3)  and  (B4) 
respectively.  The  expressions  for  Z,,  COS  <I>,,  and  SIN  O,  are 
given  by  (B6) ,  (B7) ,  and  (B8) ,  respectively.  The  dimensionless 
parameters  are  listed  in  (39)  and  a  new  dimensionless  parameter 
is  introduced: 

For  small  values  of  y,  small  angle  approximations  can  be 
employed  and  higher  order  terms  can  be  neglected.  Because  the 
extension  of  the  two  "legs"  of  the  global  bending  mechanism  is 
neglected,  A*  is  limited  to  A*  ^  0.5.  Over  this  range  of  A*,  y 
£  5  degrees  for  values  of  H*  and  L*  which  are  reasonable  for  ship 
scantlings.  Applying  small  angle  approximations  and  neglecting 
higher  order  simplifies  the  following: 


cos<Dp  =  &   +  (y)  2  ■  v/2(Y)  +  (Y)1 


2 


sin*p  =  &    -    (Y)  2  +  (Y)  2 


Z„  = 


-,  1 

/2  {  ±    -   2Y  +  2S/2Y  2  ) 


p  11 

{  Y  +  2Y2  }2 
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y/2     +      (L* 

-  2) (y)  2 

2( 

i 
v^2(Y)  2    + 

l 

y)  2 
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4.3.5.  STRAIN  HARDENING:  Strain  hardening  effects  may  be 
included  in  the  global  bending  model  in  the  same  manner  in  which 
they  were  included  for  the  local  denting  model  in  section  3.3.5. 
For  the  global  bending  model,  PB/Mo  ^n  equation  (47)  is 
represented  by  (73) . 

4.4.  COMBINED  BENDING-DENTING  MODEL:  Figures  16  and  17 
show  the  combined  local  denting  and  global  bending  load- 
deflection  curves  for  typical  parameters  for  a  40,000  DWT  USDH 
tanker.  Figure  16  is  the  load-deflection  curve  for  a  =  7C/12. 
The  curve  for  a  =  n/8  is  shown  in  Figure  17.  Both  Figure 
indicate  transitions  from  the  local  denting  mode  to  the  global 
bending  mode  where  the  two  solutions  cross. 
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5.   CRUSHING  EXPERIMENTS 

5.1.  BACKGROUND:  Four  extruded  aluminium  double  hollow 
sections  were  crushed  for  this  study.  Two  experiments  were  set- 
up to  initiate  the  local  denting  mode.  The  other  two 
experiments  were  designed  to  simulate  the  global  bending  mode. 
Appendix  C  contains  descriptions  of  the  experiments. 

The  depth-to-thickness  ratio  for  these  specimens  was 
H/t„  =  20. 3,  which  is  at  the  lower  end  of  the  thin-wall  member 
spectrum.    Also,  the  material  exhibited  significant  strain- 
hardening  properties.   Appendix  D  contains  the  specification  of 
the  aluminium  specimens  that  were  crushed. 

The  extruded  aluminium  specimens  were  provided  by  Reynolds 
Aluminium,  Holland  B.V.  This  size  of  extruded  shape  was  used  in 
order  to  maintain  a  small  scale  which  could  easily  be  handled 
with  standard  laboratory  machines.  The  specimens  were  300  mm 
long,  127  mm  wide,  and  50.8  mm  deep.  The  issue  of  structural 
discontinuity  and  inhomogeneity  due  to  welds  was  avoided  by 
using  extruded  forms. 

The  double  hollow  profile  was  chosen  in  order  to  provide 
some  restraint  against  in-plane  motion  in  the  transverse 
direction.  A  multiple  cell  section  such  as  that  shown  in  Figure 
5  would  have  been  an  ideal  choice.  The  double  hollow  form, 
however,  was  selected  due  to  availability. 
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5.2.  LOCAL  DENTING:  Figures  18  and  19  show  the  load- 
deflection  curves  for  the  two  denting  experiments.  A  photograph 
of  the  deformed  specimen  from  the  first  denting  experiment  is 
shown  in  Figure  20. 

5.2.1.  DESCRIPTION  OF  DEFORMATION:  Due  to  the  significant 
strain  hardening  effects,  it  was  difficult  to  see  the 
propagation  of  the  deformation  during  the  crushing  process.  The 
photograph  in  Figure  20  shows  the  permanent  deformation  of  the 
first  dented  specimen  with  one  of  the  outer  webs  cut  away  to 
expose  the  interior  web.  The  pattern  of  deformation  is  very 
similar  to  that  postulated  for  the  simplified  model.  The  angle 
a  between  the  flange-web  connection  and  the  first  hinge  line  is 
approximately  15  degrees.  The  shape  of  the  lower  hinge  line  is 
parabolic  rather  than  straight .  The  longitudinal  length  of  the 
deformed  region,  however,  is  approximately  twice  the  depth  of 
the  web . 

In  both  denting  experiments,  fracture  occurred  at  the 
juncture  of  the  outer  two  webs  and  the  upper  flange .  During  the 
crushing  process,  the  location  of  fracture  was  obscured  by  the 
rigid  punch.  Therefore,  the  crushing  depth  (A*)  at  which 
fracture  began  was  not  recorded. 

Both  denting  specimens  showed  initial  signs  of  fracture 
along  the  second  hinge  line.  This  appears  to  be  a  result  of  the 
fact  that  a  structure  cannot  assume  infinitely  small  radii  of 
curvature.   Because  the  local  denting  model  employs  stationary 
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hinges  and  is  formulated  for  only  a  discontinuous  deformation 
field,  this  feature  is  not  included  in  the  model. 

The  ends  of  the  specimens  for  both  denting  experiments  were 
clamped  to  the  work  surface  only  at  the  lower  flange.  A  fully 
clamped  condition  was  not  achieved.  A  small  degree  of  "pull-in" 
in  the  longitudinal  direction  was  observed  during  the  crushing 
process.  The  local  denting  model  solution  was  developed  with 
the  assumption  of  fully  clamped  ends.  The  amount  of  energy 
dissipated  by  stretching  in  the  experiments  was  probably  less 
than  that  predicted  by  the  simplified  model. 


5.2.2.  LOAD-DEFLECTION  CURVES:  The  load-deflection  curves 
labeled  "MODEL"  in  Figures  18  and  19  were  developed  from 
Equation  (A5)  with  the  parameters  listed  below.  The  curves 
shown  are  the  envelopes  of  the  curves  produced  by  (A5)  .  Each 
specimen  was  modelled  as  three  structural  units  with  the 
following  geometric  parameters: 

a  =  tc/12      H*  =  20.3 
B*  =  16.9       t*=1.0 

Although  the  two  outer  structural  units  are  channel  type 
sections  rather  than  I-beam  type  sections,  it  is  assumed  that 
the  two  outer  webs  deform  in  the  same  manner  as  the  inner  web 
during  the  crushing  process. 
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Strain  hardening  effects  were  included  as  discussed  in 
section  (3.3.3) .  The  stress-strain  curve  for  the  specimens 
which  was  provided  by  the  manufacturer  is  included  in  Appendix 
D. 

As  shown  in  Figures  18  and  19,  the  local  denting  model 
over-predicts  the  energy  absorbed  in  the  crushing  process .  It 
also  does  not  capture  the  load  reduction  which  occurs  at 
approximately  A*  =  0.10.  The  model  shows  an  almost  linear 
increase  in  load  after  the  initial  phase.  Some  of  the  load 
reduction  may  be  due  to  the  lack  of  fully  clamped  end 
conditions . 

From  the  derivations  in  Appendix  A,  it  can  be  shown  that 
the  membrane  energy  involved  in  stretching  the  flange  is  the 
most  significant  factor  in  this  region  of  the  load-deflection 
curve.  The  relative  magnitude  of  the  crushing  load  for  the 
model  and  the  general  character  of  the  curves,  however,  show 
good  agreement  with  the  experimental  results . 


5.3.  GLOBAL  BENDING  MODE:  Figures  21  and  22  show  the 
load-deflection  curves  for  the  two  bending  experiments .  Figure 
23  is  a  photograph  of  the  deformed  specimen  from  the  first 
bending  experiment . 
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5.3.1.  DESCRIPTION  OF  DEFORMATION:  Again,  the 
photograph  in  Figure  23  shows  the  permanent  deformation  viewed 
after  experiments  were  completed.  The  deformation  is  similar  to 
that  postulated  for  the  bending  model .  The  angle  a  is 
approximately  22  degrees.  The  parabolic  shapes  of  the  hinge 
lines  are  more  pronounced  in  these  specimens  than  in  those  of 
the  denting  experiments. 

In  the  first  bending  experiment,  fracture  occurred  at  the 
juncture  of  the  two  outer  webs  with  the  upper  flange .  The 
second  bending  experiment  was  halted  at  a  smaller  value  of  A*. 
The  initiation  of  fracture  is  evident  at  the  same  locations  for 
this  specimen. 

5.3.2.  LOAD-DEFLECTION  CURVES:  Figures  21  and  22  show  the 
load-deflection  curves  for  the  bending  experiments.  The  curves 
labeled  "BENDING  MODEL"  and  "DENTING  MODEL"  were  developed  from 
the  models  discussed  in  sections  3  and  4  with  strain  hardening 
effects  included  as  discussed  in  sections  3.3.5  and  4.3.5. 

The  specimens  were  modelled  as  three  structural  units  with 
the  following  geometric  parameters : 

a  =  nl    8     H*  =  20.3 
B*  =  16.9      t*  -  1.0 
L*  =  4.4 
The  global  bending  mode  was  treated  as  bending  of  a  free 
beam.   The  contributions  from  the  secondary  hinge  mechanisms 
were  neglected. 
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The  theoretical  model  for  the  bending  mode  shows  very  close 
agreement  with  the  experimental  results  both  in  the  shape  and 
the  magnitude  of  the  load-deflection  curves .  The  theoretical 
denting  model,  however,  over-predicts  the  required  force  during 
the  initial  phase  prior  to  transition  to  overall  bending. 
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6.   CONCLUSIONS  &  RECOMMENDATIONS 

The  objective  of  this  study  was  to  investigate  the  crushing 
characteristics  of  web  girders  for  a  USDH  structure.  The 
crushing  experiments  that  were  conducted  provide  more  insight 
into  the  shape  of  the  deformation  that  a  web  girder  undergoes 
when  subjected  to  a  concentrated  lateral  collapse  load.  The 
shapes  of  the  deformed  specimens  shown  in  Figures  20  and  23  are 
very  similar  to  the  photograph  in  Figure  6  which  shows  damage  to 
a  fully  clamped  plate  subjected  to  a  crushing  load. 

The  local  denting  model  over-predicted  the  amount  of  energy 
absorbed  for  the  denting  experiments.  A  model  with  a  continuous 
deformation  field  may  better  predict  the  initial  phase  for  the 
denting  solution. 

From  the  derivations  of  the  load-deflection  for  the  denting 
model,  it  can  be  seen  that  the  energy  involved  in  stretching  the 
flange  is  the  most  significant  term  contributing  to  the  total 
rate  of  internal  energy  dissipation.  The  bending  energy  and  the 
energy  involved  in  stretching  the  web  are  negligible  in 
comparison.  Denting  experiments  with  fully  clamped  end 
conditions  would  provide  better  insight  into  the  stretching 
involved  in  deforming  the  flange  plating.  Also,  experiments 
using  thinner  specimens  (greater  depth-to-thickness  ratios)  with 
less  significant  strain  hardening  characteristics  would  provide 
more  insight  into  the  initial  mode  of  deformation  and  any 
transition  to  different  modes  of  deformation. 
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The  theoretical  load-deflection  curve  for  the  global 
bending  model  closely  followed  the  experimental  results  for  a 
structural  unit  with  free-free  end  conditions.  The  shape  of  the 
load-deflection  curve  for  the  global  bending  model  is  very 
similar  to  the  shape  of  experimental  results  presented  by  de 
Oliveira  [12)  for  a  horizontally  free  circular  tube  subjected  to 
a  concentrated  lateral  crushing  load. 

The  membrane  energy  involved  in  the  extension  of  the  two 
"legs"  of  the  bending  mechanism  was  neglected  for  the 
theoretical  model.  This  may  be  a  significant  factor  for  a  fully 
clamped  structural  unit.  Crushing  experiments  for  fully  clamped 
conditions  would  be  worthwhile. 

The  models  presented  in  this  paper  were  based  on  stationary 
straight  hinge  lines.  Although  the  actual  hinge  lines  are 
parabolic  curves,  the  models  provide  reasonable  estimates  of  the 
energy  absorbed  during  the  crushing  process.  For  a  given  set  of 
geometric  parameters  (H*,  B*,  T*,  and  L*)  ,  the  formulas  developed 
for  these  models  could  be  approximated  by  simple  functions  and 
then  incorporated  in  damage  prediction  models.  For  example,  a 
linearly  increasing  function  could  be  used  to  approximate  the 
local  crushing  solution.  For  the  combined  denting-bending  mode, 
a  mean  crushing  force  could  be  employed.  These  simple  functions 
could  be  used  to  describe  the  springs  of  the  USDH  model  shown  in 
Figure  3.  The  denting  model  could  be  used  to  describe  the 
vertical  springs .  The  global  bending  model  could  be  used  to 
describe  the  diagonal  springs . 
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For  this  study,  several  assumptions  were  made  in  order  to 
keep  the  analysis  of  the  crushing  characteristics  relatively 
simple  and  mathematically  tractable.  These  assumptions  may 
require  further  investigation. 

In  particular,  the  assumption  that  there  is  no  load- 
interaction  between  bending  moments  and  membrane  forces  was  made 
to  uncouple  the  energy  dissipation  terms.  The  inclusion  of  the 
interaction  surface  for  the  different  forces  would  make  the 
problem  significantly  more  difficult,  as  pointed  out  by  de 
Oliveira  [131  . 

A  further  refinement  of  the  models  presented  in  this  paper 
could  be  made  for  the  method  of  including  strain  hardening 
effects.  Strain  hardening  was  accounted  for  by  a  relatively 
crude  manner  which  utilized  the  average  strain  over  the  entire 
structure .  The  strain  for  each  deforming  region  could  be  used 
to  adjust  the  flow  stress  for  each  region. 

Another  factor  which  requires  further  investigation  is 
fracture.  Not  only  may  fracture  affect  the  deformation  mode  and 
the  load  carrying  capacity  of  a  web  girder  unit,  it  may  also 
initiate  flooding  in  a  ship's  double  hull.  For  the  models,  an 
estimate  of  initial  fracture  could  be  made  based  using  a 
critical  strain  to  rupture  criteria  for  the  regions  of  maximum 
strain. 

For  the  crushing  experiments,  fracture  first  occurred  at 
the  upper  web  flange  connections  for  the  outer  flanges.  The 
ultimate  strain  for  the  experimental  specimens  was  on  the  order 
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of  20%.  Strain  in  this  connection  region  should  be 
investigated.  The  extruded  shapes  had  stress  concentrations  in 
this  region  due  to  the  geometry  shown  in  Figure  2£.  A  USDH  ship 
structure,  also,  will  have  stress  concentrations  due  to  welded 
joints . 
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APPENDIX  A:   LOCAL  DENTING  MODE 

In  this  section  the  development  and  the  normalization  of 
the  load-deflection  relationship  for  the  local  denting  mode  is 
detailed. 

A.l.  LOAD-DEFLECTION  RELATIONSHIP:   Equating  the  external 
rate  of  work  to  the  internal  rate  of  energy  dissipation  gives: 


»       «        »»  ♦  .  _  -  . 

Wext      =  E total      ~      EB      +      EM  flange      +      EM  web  ^JU"' 


Substituting  (37)  ,   (19)  ,   (28)  ,  and  (35)  into  (Al)  and  making 
dimensionless  yields: 

PB  £*2  2  (C3  -  C4  cos*  -  COS2*) 

M0         (£*<!  +  A*2)      cosasin*  yJ^+C^os^b-cos2^ 

+  4J2CH*   r   (cosacosp-cos2acos4>) 

{  cosP  (sin*-cos*)  +cosa  }2 

sina  (cos4>-sinfl>)  ■, 


cosa  (l+2cos*sin$) 


*  3J2  B<t<*       (costt-sintt)    }         (A2) 
(l+2cos$sin<D) 


where 


52 


cos*  =  *± 


72   (C!  +  A*) 
2 


sin*  =  Jt± 


7C*2  +  A*2 
(C*  -  A*) 


vT 


*2 


♦  2 


(A3) 


The  dimensionless  parameters  variables  are: 


A*  =  A  .   C*  =  1 


s*  =  JL  ;   t*  -  -He 


T  ' 


(A4) 


Substituting    (A3)    into    (A2)    yields: 


a   _ 


Mr 


1   +  i|r 


2\ 1/2- 


2  [iq  (l+i|r2)  -2i|f-iC2  (1+iJj)  (l+i|/2)1/2] 
v    cosa(l-i|r)^3(l+i|r2)  -2i|r+JC4  (l+i|r)  (l+i|f2)1/2; 


+    2/2  C*  H* 
1   +  i|r2 


'  [K5  (l+i|r2)  -J^(l+i|r)  (l+ijf2) 1/2]  ' 


[J^ilr2-*^  (l+lr2)  1/2+cosa] 
f  (4  C'ff'tana  +  8  B't*2)  ) 


1   +  i|rJ 


i|r    (1   +  t|/2)  1/2 


(A5) 


*i  = 

*2  = 

*3  = 

*4  = 

*5  = 

*6  = 

*7  = 


2  cos2a  +  2^2  cosacosp  +  1 
2  cos2a  +  <J2  cosacosp 
1-2 (cos2a  -  cos2p) 

2v/2  cosacosp 
2  cosacosp 
sj2  cosa 
2  cos2P  +  cos2a 


(A6) 
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i|r  = 


A! 


(A7) 


If  y  <  1,  then  y2  «  1 .   If  this  case  holds,  then  Equation 
(A5)  can  be  linearized  by  neglecting  Xj/2 .   Equation  (A5)  becomes: 


B 


Mr 


2  [    {Kx   -  K2)    -    (iq   +  2)!|r  ] 


COSO(l    -   i|i)    ^    (K2    -KJ    +    (ff4   -  2)i|r    J 

h  ^272    C*H*    [    (K5    -  K^)    -  K^ty   )  y 
cos2a  -  K^ty 

+  (  i4r    [4  C#*  tana   +  8  B*t*2]    ) 


(A8) 
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APPENDIX  B:   GLOBAL  BENDING  MODE 

B.l.  BENDING  ENERGY  PATE:   Substitution  of  (60),  (61),  and 
(62)  into  (11),  (12),  and  (13)  gives  expressions  for  [9j  ,  [62] , 
and  [03]  in  terms  of  y.   Substitution  of  these  expressions  and 
(55),  and  (56)  into  (53)  yields: 


»  ,  (  C-c^cosfc^-cos2* 


-3    ^!_l  p 


\ 


cosasin4p  yJC1+c2cos<bp-cos2<bp/ 
+     2  BW0y 


(Bl) 


where  Cx,    and  C2   are  given  by    (14);      C3,and  C4   are   given  by    (20); 
and: 


p 


2  {sin2YV^in^T~sinr2Y+cos2Y  (  —  -sin2y) )  (B2) 


2 


{sin2Y-sin22y}  2{l-sin2Y+2sin22Y+2cos2Yv/sin2Y-sin22Y) 2 


cos$p  =  -i_{  (l-sin2Y)    +   (cosY+sinY)  \/sin2Y  }  (B3) 


sinfl>p  =  -£—  {  cos2y  -    (cosY-sinY)  \/sTn2Y  }  (B4) 
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B.2.  SECONDARY  HINGE  MECHANISMS:  The  bending  energy  rate 
terms  and  the  membrane  energy  rate  terms  for  the  secondary 
mechanisms  are  equal  to  those  for  the  primary  mechanism 
evaluated  at  y/2  rather  than  y.  The  bending  energy  rate  term 
is : 


EBg=  2  ■?  M0  H  Zt 


C3-c4cos$s-cos2*, 


cosasin*s  ^C1+c2cos4s-cos2<&S; 
BW0y 


<B5) 


where 


Zs  = 


* 


2  {  sinyVsiny  -  sin2y   +  cosy(—   -  siny)   } 
{siny-sin2y  }2  {  l-siny+2sin2y+2cosy\/sinY-sin2Y  }2 


(B6) 


cos$s  =  2—  {   (1  -  siny)  +  (siny  +  sin2y)  2 


(B7) 


sin$s  =  i—  {  cosy  -   (siny  -  sin2y)  2  } 


(B8) 
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The  membrane  energy  rate  term  for  the  secondary  mechanisms 


is: 


4.  =   V2   i9N0H2 


cosacosp  -  cos2acos4s 
{  cosP  (sin*     -  cos*  )    +  coso  }2 


(B9) 


B.3.  LOAD DEFLECTION  RELATIONSHIP:    The  total  internal 

rate  of  energy  dissipation  consists  of  the  energy  rate  terms  for 
the  one  primary  mechanism  and  the  energy  rate  terms  for  the  two 
secondary  mechanisms.  Combining  these  terms,  equating  the  sum 
to  the  external  rate  of  work,  and  making  dimensionless  yields: 


s_ 


Mn 


4Z, 


C3-C4cos*p-cos2$J5 


2y/2H*  (cosacosp-cos2acos*s) 
cosasinG^q  +  C^cosGp-cos2*^       (cosP  <sin<&p-costfp)  +cosa) 2 


BZ, 


C3-C4cos*s-cos2$, 


2y/2H*  (cosacosp-cos2acos*s) 
cosasin*sV/C1  +  C2cos*s-cos2*s      (cosP  (sin*s-cos*s)  +cosa)2 


4B*t*3 
ZjH* 


(BIO) 
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where 


i 
A*  =  72  (Y)  2 


'^1-iIy  <b11> 


The  dimensionless  parameters  are  listed  in  (39)  and  a  new 
dimensionless  parameter  is  introduced: 

L*  =  -  (B12) 

H 
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APPENDIX  C:   EXPERIMENT  SET-UP 

C.l.  DENTING  EXPERIMENTS:  Figure  24  shows  a  side-view  and 
an  end-view  of  the  set-up  for  the  denting  experiments.  The 
specimens  were  clamped  to  the  work  surface  only  along  the  lower 
flange.  The  work  surface  (the  platen)  was  raised  in  a 
displacement  control  mode  and  the  reaction  force  PB  was  measured 
as  a  function  of  displacement.  A  solid  steel  cylinder  25  mm 
diameter  was  used  as  the  rigid  punch.  The  solid  cylinder  was 
used  rather  than  a  knife-edge  indenter  to  delay  initiation  of 
fracture . 

C.2.  BENDING  EXPERIMENTS:  Figure  25  shows  a  side-view  and 
an  end-view  of  the  set-up  for  the  bending  experiments.  A  three- 
point  bending  set-up  was  used.  The  specimens  were  not  clamped. 
Instead,  bending  of  a  free  beam  was  modelled.  The  specimens 
were  placed  on  two  steel  rods  of  25  mm  diameter  that  were  spaced 
225  mm  apart.  A  third  steel  rod  was  connected  to  the  load-cell 
above  the  specimens.  The  base  was  raised  in  a  displacement 
control  mode  and  the  reaction  force  PB  was  measured  as  a  function 
of  displacement . 
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APPENDIX  D:   EXPERIMENT  SPECIMEN  SPECIFICATIONS 

D.l.  DIMENSIONS:  Figure  2  6  shows  a  cross-section  and  the 
dimensions  of  the  prismatic  specimens  used  for  the  crushing 
experiments . 

D.2.  MATERIAL  PROPERTIES:  The  specimens  used  for  the 
crushing  experiments  were  extruded  60  61  aluminium  which  were 
supplied  by  Reynolds  Aluminium  Holland,  B.V.  The  material 
properties  for  the  specimens  are  listed  in  Table  D-l,  below. 

TABLE  D-l:   MATERIAL  PROPERTIES 


YOUNG'S  MODULUS 

69   GPa 

YIELD  STRESS 

114   MPa 

ULTIMATE  STRESS 

220   MPa 

ULTIMATE  STRAIN 

0.206 

BRINELL  HARDNESS 

69 

The  material  has  significant  strain  hardening  properties. 
In  order  to  incorporate  this  characteristic  in  the  theoretical 
models,  the  stress-strain  curve  was  modelled  as: 


a     =    a 


<      € 


ult 


V  euic 


(Dl) 


with  n  =   0.165;    a0  =   220   MPa;    and   eult   =    0.206 
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The  points  plotted  in  Figure  27  were  obtained  from  the 
tensile  stress-strain  curve  provided  by  the  supplier.  The  curve 
shown  in  Figure  27  is  the  power-law  model  which  was  used  for 
the  theoretical  model  solutions. 
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FIGURES 


Figure   1:      The   USDH   Hull   Form 
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Figure  2:   Transverse  Frame  Crushing  Model 


Figure  3 :   Spring  Model  of  USDH 
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Figure  4:   USDH  Section 


Figure  5:   Structural  Unit 


64 


Figure  6:   Deformation  of  Clamped  Plate 
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Figure  7:   Local  Denting  Model 


Figure  8:   Deforming  Web  Girder 


66 


Figure  9:   Deformed  Flange 
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Figure  10:   Stretching  in  Local  Denting  Mode 
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Figure    11:       Local   Denting  Mode    for   a  =   7t/12 
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LOCAL  DENTING  MODEL  (ALPHA  =  PI/8) 
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Figure  12:   Local  Denting   Mode  for  a  =  rt/8 
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Figure  13:   Three  Hinge  Mechanism 


Figure  14:   Hinge  Point  Mechanism 
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Figure  15:   Stretching  in  Hinge  Point  Mechanism 
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Figure   16:      Global   Bending  Mode    for  a  -   tc/12 
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Figure  17:   Global  Bending  Mode  for  a  -  rc/8 
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Figure  18:   First  Denting  Experiment 


75 


DENTING    #2 

CH*=20.3j    B*=16.9j    T*=1j    ALPHA   =    PI/ 12} 


200     - 


3 
2 


fi 


MXCL 


Figure  19:   Second  Denting  Experiment 
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Figure  20:   Photograph  of  Dented  Specimen 
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Figure  21:   First  Bending  Experiment 
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Figure  22 :   Second  Bending  Experiment 
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Figure  23 :   Photograph  of  Bent  Specimen 
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Figure  24:   Set-up  for  Denting  Experiments 


Figure  25:   Set-up  for  Bending  Experiments 


81 


1% 


uoT^oes-ssoao   uauiTosds      -SZ   •atnBfj 


<S> 


(X8)5'9  y- 


V/W/V/y  ■  zzzzzzzz 


^ 


K 


/ 


i&ss^zzz 


J5 


(X8i  f  a 


2ZZZZZZZZZZZZZZZZZ 


& 


Q'2 


TT5T 


STRESS -STRAIN  CURVE  FOR  ALUMINIUM  SPECIMENS 

n  -  0.165 


1 

■          > 

j» — - 

O.f 

w 

in 

HI 

a 

to 

UJ 

n.B 

K 

I 

I- 

-I 

D 

0.4 

•^ 

\n 

in 

S 

K 

</) 

O.g 
0  » 

I L_ 

■                        i                        i 

0  4  0  1 

STRAIN/    ULTIMATE    STRAIN 
«    actus  I   vo I ues 


Figure  27:   Stress-Strain  Curve  for  Specimens 
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